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Abstract. An introduction and survey is given of some recent 
work on the infinitesimal dynamics of crystal frameworks, that is, 
of translationally periodic discrete bond- node structures in M'', for 
d = 2, 3, . . . . We discuss the rigidity matrix, a fundamental object 
from finite bar-joint framework theory, rigidity operators, matrix- 
function representations and low energy phonons. These phonons 
in material crystals, such as quartz and zeolites, are known as 
rigid unit modes, or RUMs, and are associated with the relative 
motions of rigid units, such as Si04 tetrahedra in the tetrahedral 
polyhedral bond-node model for quartz. We also introduce semi- 
infinite crystal frameworks, bi-crystal frameworks and associated 
multi-variable Toeplitz operators. 



1. Introduction 

A survey is given of some recent work on the infinitesimal dynamics 
of crystal frameworks, by whicfi we mean translationally periodic dis- 
crete bar-joint frameworks in W^. This includes a discussion of rigidity 
operators, matrix symbol function representations and the connections 
with models for low energy phonon modes in various material crystals. 
These modes are also known as rigid unit modes, or RUMs, reflect- 
ing their origin in the relative motion of rigid units in the crystalline 
structure. I also introduce briefly the contexts of semi-infinite crystal 
frameworks and bicrystal frameworks and indicate how their rigidity 
operators involve multivariable Toeplitz operators whose symbol func- 
tions are matrices over multi-variable trigonometric polynomials on the 
d-torus. 

The topic of infinite bar-joint frameworks, whether periodic or not, 
can be pursued as a purely mathematical endeavour and many aspects 
of deformability and rigidity remain to be understood. The main per- 
spectives below and related issues are developed in Owen and Power 
122], [22], [25] and Power [2S], [2Z]. 

Translationally periodic bond/node bar-joint frameworks or networks 
are ubiquitous in mathematics (periodic tilings for example), solid state 
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physics (crystal lattices, graphene), solid state chemistry (zeolites) and 
material science (microporous metal organic frameworks). So there is 
no lack of interesting examples. I shall illustrate a number of concepts 
with three examples derived from tilings seen in Seville at the Alcazaar 
and the Cathedral. 

2. Models for material crystals and low energy phonons. 

We begin by outlining one particualr motivation from material sci- 
ence. A crystal framework C in can serve as a mathematical model 
for the essential geometry of the disposition of atoms and bonds in a 
material crystal M.. In the model of interest to us the vertices cor- 
respond to certain atoms while the edges correspond in some way to 
strong bonds. Also the identification of strongly bonded "units" in M. 
imply a polyhedral net structure and it is this that gives the relevant 
abstract framework C. A fundamental example of this kind is quartz, 
Si02, in which each silicon atom lies at the centre of a strongly bonded 
Si04 unit, which in turn may be modeled as a tetrahedron with an 
oxygen atom at each vertex. In this way the material crystal quartz 
provides a mathematical crystal framework of pairwise connected tetra- 
hedra with a particular connectedness and geometry. 

Material scientists are interested in the manifestation and explana- 
tion of various forms of low energy motion and oscillation in materials. 
Of particular interest are the rigid unit modes in aluminosilicate crys- 
tals and zeolites, where quite complicated tetrahedral net models are 
relevant. These low energy (long wavelength) phonon modes are ob- 
served in neutron scattering experiments and have been shown to cor- 
relate closely with the modes observed in computational simulations. 
There is now a considerable body of literature tabulating the (reduced) 
wave vectors of RUMs of various crystals as subsets of the unit cube 
(Brillouin zone) and it has become evident that the primary determi- 
nant is the geometric structure of the abstract frameworks C. See, for 
example. Dove et al [8], Hammond et al [12], [13], Giddy et al [10] and 
Swainson and Dove [32] . Particularly intriguing is the simulation study 
in Dove et al [8] which gives a range of pictures of the RUM spectrum 
and multiplicites for various idealized crystal types. 

In the experiments and in the simulations the background math- 
ematical model is classical lattice dynamics and rigid unit modes are 
observed where the phonon dispersion curves indicate vanishing energy. 
However, one can also identify such limiting cases through a direct lin- 
ear approach as we outline below and from this it follows that these 
sets (at least in simulations) may be viewed as real algebraic varieties. 
See Theorem 15.41 below, [25] , [26] and Wegner [31] . It is convenient to 
define the RUM dimension to be the dimension of this algebraic variety. 
(See Section 5.) In 3D it takes the values 0, 1, 2, 3. 
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3. An ILLUSTRATIVE EXAMPLE 

The following simple example will serve well to illustrate the notation 
scheme for general crystal frameworks in d dimensions that we adopt. 
The example is also of interest in its own right, as we see later. 

Figure 1 indicates a translationally periodic bar-joint framework 
C = {G,p) determined by a sequence p = (pk) in R^. The frame- 
work edges [pi,pj], associated with the edges of the underlying graph 
G, are viewed as inextensible bars connected at the framework vertices 
Pk but otherwise unconstrained. 




Figure 1. An infinite bar-joint framework. 

Let the scaling be such that 

Pi = Plfi = (0, 0), P2= P2,0 = (0, 4), P3 = p3,o = (1, 3) 

are three framework vertices of a triangular subframework. Write their 
translates as 

PK,k = TkP^fl, for K G {1,2,3}, A; e Z, 

where is the isometry : (x, y) — >■ {x+Ak, y) od M^. The translation 
group T = {Tfc : k E Z,} is also used to define a natural periodic 
labelling of the framework edges: 

ei = ei,o = [Pl,P2], 62 = 62,0 = [P2,P3], 63 = 63,0 = [Pl,P3,], 
64 = 64,0 = [Pl,0,Pl,l], 65 = 65,0 = [p3fi,P2,l] 

and 

Cj- fc = TkCj, , for j e {I, 2, 3, 4, 5}, k e Z. 
Thus, the pair of finite sets 

Fv = {Pl,P2,P3}, Fe = {ei, . . . ,65} 

have disjoint translates under T and the set of all such translates de- 
termines C. 

In an exactly similar way a translationally periodic bar-joint frame- 
work C in R'^ is determined by a triple [Fy, F^, T) where we refer to 
the finite set pair Ai = {F^, F^) as a motif for C. Of particular interest 
for applications are the cases d = 2, 3 in which T = {T^ : k G Z'^} and 
T has full rank. For d = 3 " full rank" means that the so-called period 
vectors 

0-1 = F(i,o,o)(0), a2 = T(o,i,o)(0), = T(o,o,i)(0) 
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are linearly independent, in which case the framework vertices, if they 
are distinct, form a discrete set in M^. We call such discrete bar-joint 
frameworks crystal frameworks. 

We now introduce a key dynamical ingredient, namely the notion of 
an infinitesimal flex. This definition is the same as that for finite bar- 
joint frameworks being a specification of velocity vectors at the nodes 
which, to first order, do not change edge lengths. 

Definition 3.1. Let C be crystal framework, with framework vertices 
p^.fc as above. An infinitesimal flex of C is a set of vectors (velocity 
vectors) such that for each edge e = [pK,k,PT,i] 

{PK,k - Pt,1, Un,k - Ur,l) = 0. 

The linear equations required for an infinitesimal flex u = (u^^k) 
translate to a single equation R{C)u = where R{C) is the so called 
rigidity matrix for the framework and where m is a vector in the direct 
product vector space "H^ = YIk k ^"^5 regarded as a composite vector 
of instantaneous velocities. The rigidity matrix is sparse with rows 
labelled by edges and columns labelled by the Euclidean coordinate 
labels (k, k, a) of the framework vertices, with a G {1, . . . ,d}; the row 
for framework edge e = [pi,Pj] has the entry (p^./c — Pt,i)u for column 
(k, k, cr), and has the negative of this entry for column (r, /, a). Thus 
for (i = 3 row e appears as 

[0 ... t^e ... - . . . 0] 

where the vector Ve = PK,k—Pr,i (resp —Ve) is distributed in the columns 
for {K,k,a) (resp. (r, Z,(j)) with a G {x,y,z}. 

From various viewpoints, such as phase-periodic velocity vectors on 
the one hand, or square-summable velocity vectors on the other hand, 
with the introduction of complex scalars and functional representations 
of vector spaces, the rigidity matrix R{C) leads to a matrix-valued 
function $(2) with |Fe| rows and d\Fe\ columns. The entries are scalar- 
valued functions on the (i-torus of points z = {zi,. . .Zd) in with 
= 1. 

We define this matrix function below in Definition 15.31 and one can 
check that the strip framework of Figure 1 has associated matrix func- 
tion 
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3.1. Examples from Seville. The next three frameworks are based 
on some simple two-dimensional tessellations that are suggested by 
tilings found in Seville cathedral (Figures 2 and 3) and in the Alcazaar 
in Seville (Figure 4). All three are in Maxwell counting equilibrium 
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in the sense that the average coordination number (average vertex de- 
gree) is 2, matching the degree of freedom of each vertex, while each 
subframework is not " overconstrained" , in the sense that the number 
of edges does not exceed twice the number of vertices. 

A motif for the framework Cgevi is shown in Figure 3 together with 
the period vectors (dotted). The motif edges consist of a square of 
edges together with two vertical edges whose (equal) lengths fix the 
geometry up to a global scaling. (The "rigid units" of this framework 
are the single vertex subframe works.) 

Note that there is an evident (edge-length preserving) continuous 
flex, or deformation, p{t) = {pK,k{t)) of Cgcvi which is associated with 
an expansion in the x direction and a matching contraction in the y 
direction. We remark that in the case of the geometry with all edge 
lengths equal this deformation passes through the framework composed 
of congruent rhombs which is reciprocal (in the lattice sense [6]) to the 
well-known kagome framework, indicated in Section 5.1. From the first 
instant of the deformation, so to speak, one obtains an infinitesimal flex 
u = p'{0) which (unlike infinitesimal translation flexes) is unbounded. 
Less evident are various nontrivial bounded infinitesimal flexes, but we 
see below that there are plenty of these and that the RUM dimension 
is 1. 




Figure 2. The crystal framework C, 



The framework Csev2 in Figure 3 has triangular rigid units and a local 
(finitely supported) infinitesimal flex is indicated, with four nonzero 
velocity vector components. It is a general principle, as we note further 
below, that such a local phenomenon makes the framework maximally 
flexible from a RUM point of view. 

The planar graph, or "topology" of this framework is rather interest- 
ing, being a network of 4-rings of triangles connected " square- wise" , so 
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Figure 3. Csev2! with a local infinitesimal flex. 



that the "holes" are 8-cycles and 4-cycles. One can make similar con- 
structions with equilateral triangles with this topology although now 
the local infinitesimal flex is lost. We remark that periodic networks of 
pairwise corner-joined congruent equilateral triangles provide the 2D 
variants of the tetrahedral nets associated with zeolites. (That the hole 
cycles of 2D zeolites can be arbitrarily large follows from the substitu- 
tion move which replaces each rigid unit triangle by a 3-ring of smaller 
triangles with edge length halved.) 




Figure 4. Part of Csevg, a homeomorph of C12. 



The framework Cgevs is derived from a tiling in the Alcazaar in Seville. 
A moment's thought reveals it to have the same topology (underlying 
graph) as that for the basic square grid framework Cj2. The infinitesi- 
mal flexibility is less evident than is the case for Cjp. but from the RUM 
viewpoint they turn out to be equally flexible with RUM dimension 1. 
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4. Bar-joint frameworks - a very brief overview. 

4.1. Watt, Peaucellier, Cauchy, Euler, Kempe, Maxwell, Laman 

Informally, a "linkage" is a bar-joint framework with one degree of es- 
sential flexibility. In 1784 James Watt designed a bar-joint linkage 
which transformed circular motion into approximate linear motion. 
This was rather important for steam engine transmission. The mecha- 
nism was approximate and was superseded by Peaucellier's exact linear 
motion linkage eighty years later. In 1876 Kempe [19] solved the gen- 
eral inverse problem by showing that any finite algebraic curve can be 
simulated by a linkage. The rigidity of geometric frameworks was also 
of interest to Euler and to Cauchy who in particular were concerned 
with the rigidity of polyhedra with hinged faces; a beautiful classical 
result is the infinitesimal rigidity of all convex triangle- faced polyhedra. 

James Clerk Maxwell initiated combinatorial aspects with the obser- 
vation that a graph G = (V, E) with a minimally rigid generic frame- 
work realisation in the plane must satisfy the simple counting rule 
2\V\ — \E\ = 3 together with the inequalities 2\V'\ — \E'\ > 3 for all 
subgraphs. The number 3 represents the number of independent global 
infinitesimal motions for the plane. In 1970 Laman [2U] obtained the 
fundamental result that Maxwell's conditions are sufficient for generic 
rigidity and this result also anticipated the advent of matroid theory 
in rigidity theory. While corresponding counting rules are necessary in 
three dimensions they fail to be sufficient and no necessary and suf- 
ficient combinatorial conditions for generic rigidity are known ! For 
further information see [I], [2], [11]. 

4.2. Some Recent work. Laman's theorem concerns generic frame- 
works with a particular graph. One can expect that special frameworks 
with global symmetries may have more flexibility and this is a topic of 
current interest. See for example, Connelley et al [S], Owen and Power 
[23] and Shulze [22], [3D]- Understanding constraint systems of geomet- 
ric objects with symmetry present is also of significance for algorithms 
for CAD software [23] . 

Laman's theorem is also concerned with finite frameworks. A nat- 
ural generalisation, also of significance for applications, are periodic 
frameworks in the sense of the crystal frameworks above. See also 
Whiteley [33], Borcea and Streinu [1], Malestein and Theran [2T], and 
Ross, Schulze and Whiteley [2B] . 

The theory of general infinite bar-joint frameworks, from the point of 
view of rigidity, is a novel topic and perhaps a rather curious one. We 
point out in Owen and Power [21], [22] that it is possible to generalise 
Kempe's theorem to the effect that any continuous curve (ie. continu- 
ous image of [0, 1] in M^) may be simulated by an infinite linkage. (In 
[21] this is achieved with three vertices of infinite degree but in fact 
infinite degree vertices are not necessary.) 
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In material science the microporous flexing materials known as ze- 
olites are on the one hand important for industrial applications, as 
filters, and on the other hand present diverse tetrahedral rigid unit 
frameworks. The degree of continuous flexibility of such idealised zeo- 
lites is investigated in Kapko et al [18]. 

5. The RUM spectrum and RUM dimension of C 

Let C = {Fy, Fe, T) be a crystal framework in R"^ and let ICy be the 
vector space ^ C^l^^l consisting of infinitesimal velocity vectors. Let 
T*^ be the (i-torus of points u = (wi, . . . , Ud) and for /c G Z*^ write 
for the unimodular complex numbers . . . w^''. 

Definition 5.1. ([25], [26].) (a) A velocity vector u in ]Cy is periodic- 
modulo-phase for the (multi-)phase factor cu G T*^ if there exists a 
vector u = in C'-^"' such that 



Also ICy denotes the associated vector subspace of such vectors. 

(b) A periodic-modulo-phase infinitesimal flex (or wave flex) is a 
vector u in /C^ which is an infinitesimal flex for C. 

(c) The rigid unit mode spectrum, or RUM spectrum, of C (with 
specified translation group T) is the set Q{C) of phases for which there 
exists a nonzero periodic-modulo-phase infinitesimal flex. 

To each multiphase u there exists a unique wave vector k = (ki, . . . , k^) 
such that Uj = e^'^* , 1 < J < (i- Ignoring bond constraints for the 
moment, recall that the framework points of C undergo harmonic mo- 
tion with wave vector k when the vertex positions at time t satisfy 
equations of the form 



where a/27r is the frequency of the oscillation and where k ■ /c is the 
inner product ki ■ ki + ■ ■ ■ + kd ■ ka- Such pure motions appear as 
basic solutions in lattice dynamics, under harmonic approximation, 
with general solutions obtained by linear superposition. (See Dove [H] 
for example.) 

The following theorem from [26] provides an explanation for the con- 
nection between low energy oscillation modes alluded to in Section 2 
and infinitesimal wave flexes. 

Theorem 5.2. LetC be a crystal framework, with specified periodicity, 
and let k be a wave vector with point w G T*^. Then the following 
assertions are equivalent. 

(i) {uj^Uf^)f^^k is a nonzero periodic-modulo-phase infinitesimal (com- 
plex) flex for C. 

(a) For the vertex wave motion 




KG Fy,ke U^. 



PK,k{t) = p^^k + exp(27rz k ■ k) exp{iat)v, 



+ ■UKexp(27ri k ■ k) exp(2at). 
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and a given time interval, t G [0,T], the bond length changes 

for the edges e tend to zero uniformly, in both t and e, as the wavelength 
21: j a tends to infinity. 

Next we define the matrix function $c of a crystal framework C with 
given period vectors. For the multi- index k = {ki, . . . , kd) write for 
the usual monomial function on T*^. 

Definition 5.3. Let C be a crystal framework in R'^, with motif sets 

= K,o :!</«< \F,\}, = {ei : 1 < z < 

and for each edge e = [pK,k, Pt,i] in let be the edge vector p^.fc —Pt,i- 
The matrix- valued function $c(^) has rows labelled by the edges e G -Fg 
and has d\F^\ columns labelled by n and the coordinate index a G 
{1, . . . , d}. If K ^ T for edge e in Fe then 

($c(2;))e,(K,a) = {Ve)aZ'', 
($c(2))e,(r,a) = -{Ve)aZ^, 

while for each reflexive edge [pK,o,PK,5{e)] 

with the remaining entries in each row equal to zero. 

The next theorem gives one connection between ^c{z) and the in- 
finitesimal flex properties of C. Here we view the rigidity matrix R{C) 
as a linear transformation from the product vector space K.^ = Yin k '^'^ 
to the edge vector space /Ce = Hedges ^ = Ile.fe 

Theorem 5.4. The restriction of the rigidity matrix R{C) to the finite 
dimensional vector space IQ^ has representing matrix $c(^) with respect 
to natural vector space bases. 

Proof. Let m be a velocity vector in /C^ determined by u G C^l^-I as 
above. Let e in Fg be an edge of the form [pt^^k,PT,k+&\- Let (■, ■) denote 
the bilinear form on C^. The (e, A;)*'* entry of R{C)u can be written as 

{R{C)u)e,k = {Ve,UK,k) - {Ve,Ur,k+s) 
= UJ^{{Ve,U^) + (-U;V,Mr)), 

and one can check that this agrees with u;^'($c(^)'^^)e) both in the case 
K ^ T and when k, = t. □ 

We can now identify the R UM spectrum of C as the algebraic variety 
in T*^ given by 

n{C) = {z : rank$c(2) < |i^e|}. 
This set does depend on the choice of translation group. One could 
define the primitive RUM spectrum to correspond to the translation 
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group for a primitive unit cell, and this is then well-defined, up to 
coordinate permutations. If one doubles the period vectors, and hence 
the unit cell, then the new spectrum is obtained simply as the range of 
the old spectrum under the doubling map (wi, 1^2,^3 ) {wl,W2,wl). 

While we have given the multiphase form of the RUM spectrum the 
convention in material science is to indicate such a spectrum (in three 
dimensions) as the set of (reduced) wave vectors k in the unit cube 
[0, 1)'^. For calculations of RUM spectrum by different methods see 
Dove et al |8] (simulation calculations), Wegner |3l] (computer algebra 
calculations), Owen and Power and Power [SB] (direct calculations). 

The algebraic variety perspective of RUMs appears to be new and 
opens the way for new methods and terminology for understanding the 
curious curved surfaces in [5]. For example, it is natural to define the 
RUM dimension of C to be the topological dimension of Q{C) as a 
real algebraic variety. By the comments above on unit cell doubling it 
follows that this quantity is independent of the translation group. 

Tetrahedral net frameworks in 3 dimensions, with pairwise vertex 
connection, satisfy Maxwell counting equilibrium and in the periodic 
case, with no penetrating tetrahedra, are sometimes referred to as 
hypothetical zeolites [S]. (In material crystalline zeolites the rigidly 
bonded Si04 units make up such a bond-node framework.) Even in 
this case all possibilities occur for the RUM dimension, namely 0, 1, 2, 3, 
and this depends, roughly speaking, on the degree of symmetry of the 
framework. In particular as we note below the framework for the cubic 
form of sodalite indicated below has full RUM spectrum, corresponding 
to dimension 3. (This so-called order property of sodalite was first 
observed experimentally. See [12] •) 

For crystal frameworks in Maxwell counting equilibrium the matrix 
function is square and the RUM spectrum is revealed, in theory at 
least, as the intersection of the zero set of the mult i- variable polynomial 
det<I>c(-2^) with the (i-torus T'^. In fact, after fixing a monomial order 
on the d indeterminates zi, . . . ,Zd one may formally define the crystal 
polynomial pc{z), associated with C. (See [26].) This is given as the 
product z'^ det($c(-2^)) where the monomial exponent 7 is chosen so that 

(i) Pc{z) is a linear combination of nonnegative power monomials. 



(ii) pc{z) has minimum total degree, and 

(iii) pc{z) has leading monomial with coefficient 1. 

The RUM spectrum certainly has symmetry refiecting the crystallo- 
graphic group symmetries of the crystal framework. Even so the point 
group may be trivial and the following abstract inverse problem (a 
Kempe theorem for RUMS ?) may well have an affirmative answer. 
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Problem. Let q{z, w) be a polynomial with real coefficients with g(l, 1) = 
0. Is there a crystal framework with crystal polynomial p{z, w) whose 
zero set on the 2-torus is the same as that for q{z,w) ? 

5.1. Examples. The tiling-derived framework of Figure 3 has vanish- 
ing crystal polynomial. Indeed, this can be predicted from the existence 
of a local infinitesimal flex. Such a flex allows the construction of infin- 
itesimal phase-periodic flexes for all phases and so the zero set of the 
polynomial includes the entire 2-torus. 

The tiling-derived framework of Figure 2 is also in Maxwell counting 
equilibrium, its symbol function is 6 x 6 and one can show by direct 
calculation that if we write the indeterminates in this case as z, w, then 
the crystal polynomial is 

{w - l){z'^ + z{l + w) +w) 

The kagome framework is the framework Ckag formed by pairwise cor- 
ner connected equilateral triangles in regular hexagonal arrangement. 
Its symbol function is a 6 x 6 matrix and, if we write the indeterminates 
in this case as z,w, then the crystal polynomial is 

{z — l){w — l){z — w). 

There is a natural 3D variant of the kagome lattice known as the 
kagome net. The corresponding crystal framework Cknet has period 
vectors formed by three edges of a parallelapiped at pairwise angles of 
7r/3, and each parallelapiped contains two tetrahedral rigid units such 
that the three planar slices of Cknet for each pair of period vectors, is a 
copy of Ckag. The crystal polynomial takes the form 

p{z, w, u) = {z — l){w — l){u — l){z — w){w — u){z — u). 

The factorisations in these examples makes evident the nature of the 
RUM spectrum as a union of lines and a union of surfaces, respectively. 
In fact the individual factors can be predicted in terms of the identifica- 
tion of infinitesimal flexes that are supported within a linear band (for 
the 2D case) or a linear tube (in the 3D case). See [26]. For consider- 
ably more complicated polynomials with nonlinear "exotic" spectrum 
see Wegner [31] and Power [26] . 

Figure 4 shows a 4-ring of tetrahedra, three copies of which, placed 
on three adjacent sides of an imaginary cube, provide the edges and 
vertices for a motif (F^,, Ff.) for the framework Csod for the cubic form of 
sodalite. (From a mathematical perspective, this structure is arguably 
the most elegant of the naturally occurring zeolite framework types 

A full set of eight 4-rings forms a so-called sodalite cage. With the 
24 outer vertices of this cage fixed there is nevertheless an infinitesimal 
flex of the structure and so a local infinitesimal flex of Csod- This is 
in analogy with the framework Csev2- It follows that the determinant 
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of the symbol function (a 72 x 72 sparse function matrix) vanishes 
identically and that the sodalite framework Csod has RUM dimension 
3. 




6. Flexes with decay and Toeplitz operators. 

Let C be a crystal framework with an implicit choice of translational 
periodicity. Write /C^ and K.^ for the Hilbert spaces of square-summable 
sequences in K.^ and /Cg. Then R{C) determines a bounded Hilbert 
space operator from /C^ to /Cg. The natural Hilbert space basis for 
Kl, associated with the given periodicity, may be denoted {^K,a,k}, 
where a ranges from 1 to d. Similarly, the basis for K,^ is {f]e,k}^ with 
ee Fe,ke If. 

Regarding such square-summable sequences as the Fourier series of 
square integrable vector- valued functions one obtains unitary operators 

and Ue'. JCl ^ L'^{T'^) (g) Cl-^<=l. In the next theorem, from the uni- 
tary equivalence referred to is two-sided and the equivalence in ques- 
tion is the operator identity U*R{C)Uv = M$^. That U*R{C)Uv has 
the form of a multiplication operator follows from standard opera- 
tor theory, since this operator intertwines the canonical shift operators. 
Borrowing standard operator terminology we refer to $0 as the symbol 
function of C. 

Theorem 6.1. The infinite rigidity matrix R{C) of the crystal frame- 
work C in M"^ determines a Hilbert space operator which is unitarily 
equivalent to the multiplication operator 

where is the matrix function for C. 



The following corollary follows from elementary operator theory. 
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Corollary 6.2. A crystal framework has a square- summahle infinites- 
imal flex if and only if its symbol function has reduced column rank on 
a set of positive measure. 

It is natural to ask whether crystal frameworks possess infinitesimal 
flexes which decay to zero at infinity. Note that if an infinitesimal linear 
subframework has an infinitesimal fiex with such decay then the fiex ve- 
locities must be orthogonal to the direction of this subframework. (For 
otherwise there must be identical nozero velocity components in that 
direction on all the subframework points.) For this reason it follows 
that the kagome framework and similar "linear" frameworks have no 
asymptotically vanishing fiexes and in particular, no square-summable 
fiexes. In fact one can exploit the matrix function formalism to obtain 
the following much more general fact. 

Theorem 6.3. [2S] The following are equivalent for a crystal frame- 
work C with Maxwell counting equilibrium. 

(i) C has a nonzero internal ("finitely-supported" , "local") infinites- 
imal flex 

(ii) C has a nonzero summahle infinitesimal flex. 

(Hi) C has a nonzero square-summable infinitesimal flex. 

6.1. Semi-infinite and bi-crystal frameworks. We may define a 
semi-infinite crystal framework D as a subframework of a crystal frame- 
work C with an exposed face. More formally V is supported by the 
framework vertices that lie in a half-space which is invariant under a 
subsemigroup of an underlying translation group for C. In the case 
of planar frameworks this may be specified in the form of a triple 
(F^, Fe, 7+) where (F^,, Fe) is an appropriate motif and 7+ is a sub- 
semigroup of T isomorphic to one of Z_|_ x Z, Z_ x Z, Z x Z_|_, Z x Z_. 
It is not hard to verify that the rigidity operators of semi-infinite crys- 
tal frameworks can be identified with various Toeplitz operators de- 
rived from by compression to Hardy space Hilbert spaces, such as 
H'^{T) ® L2(T) in the case of Z+ x Z. 

We remark that semi-infinite frameworks have rigidity matrices that 
feature as block submatrices of the rigidity matrices of hi- crystal frame- 
works. By this we mean (for example) a framework obtained in three 
dimensions by identifying two semi-infinite frameworks at their com- 
mon surface of vertices, when this is possible. It seems that Toeplitz 
operators could provide a useful formalism for their analysisEI 

For semi-infinite crystal frameworks the equivalence of the previous 
theorem no longer holds and we illustrate this with a simple variant 
of the strip framework in Figure 1 whose matricial symbol function 
$(z) is as given in Section 2. Consider the submatrix function $o(^) 



^ Added October 2011. It appears, from [3T] for example, that there is extensive 
interest in material bicrystals. 
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obtained on removing the first two columns, corresponding to the "sup- 
porting" framework vertices, and removing the row corresponding to 
the "base" edges. The degeneracies of this matrix function correspond 
to the phases of periodic-modulo-phase infinitesimal flexes which do 
not deflect the "supporting" vertices. We have 

4 " 

-111-1 
13' 

3^ z —3 —1 

The determinant, 16(2 — 3^), does not vanish on T and so the "base- 
rooted" framework is infinitesimally rigid from the point of view of 
phase-periodic infinitesimal flexes. Since the determinant certainly 
does not vanish on a set of positive measure on T, there are no square- 
summable infinitesimal flexes which fix the baseline vertices of the 
framework. On the other hand there is an unbounded infinitesimal 
flex, corresponding to a two-way infinite geometric series for z = 2/3 
and this unbounded flex reflects the concatenated lever structure of the 
framework. 

This analysis also applies to the strip framework in Figure 6. One can 
readily check that up to scalar multiplication there is a unique proper 
(unbounded) infinitesimal flex. (Note incidentally, that this flex does 
not extend to a continuous flex. Put another way, each finite strip 
subframework here is continuously flexible but the complete two-way 
infinite framework is not. This kind of phenomenon is referred to as 
vanishing flexibility in [25] and can occur in more subtle ways.) 







,/ 


' 


/ 






• 







Figure 6. An unbounded infinitesimal flex. 



There are two natural semi-infinite frameworks C+, C_ associated 
with the strip framework of Figure 6, namely the right strip and the 
left strip. Each has a infinite triangulated rigid base framework which 
supports linked triangles. The former has a square-summable flex while 
the latter does not. This fact is evident by elementary direct analysis 
and in fact can be viewed as a reflection of the contrasting nature of 
analytic and coanalytic Toeplitz operators. We expect such operator 
theory to play a useful role in the analysis of more complex examples 
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with larger unit cells, and in the analysis of surface phonons and surface 
phenomena in semi-infinite structures. 
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